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DATA ASSIMILATION ALGORITHM FOR 3D BENARD 
CONVECTION IN POROUS MEDIA EMPLOYING ONLY 
TEMPERATURE MEASUREMENTS 

ASEEL FARHAT, EVELYN LUNASIN, AND EDRISS S. TITI 


Abstract. In this paper we propose a continuous data assimilation (down- 
scaling) algorithm for the Benard convection in porous media using only coarse 
mesh measurements of the temperature. In this algorithm, we incorporate the 
observables as a feedback (nudging) term in the evolution equation of the tem¬ 
perature. We show that under an appropriate choice of the nudging parameter 
and the size of the mesh, and under the assumption that the observed data is 
error free, the solution of the proposed algorithm converges at an exponential 
rate, asymptotically in time, to the unique exact unknown reference solution 
of the original system, associated with the observed (finite dimensional pro¬ 
jection of) temperature data. Moreover, we note that in the case where the 
observational measurements are not error free, one can estimate the error be¬ 
tween the solution of the algorithm and the exact reference solution of the 
system in terms of the error in the measurements. 
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1. Introduction 

Linking mathematical and computational models to a set of a coarse mesh ob¬ 
served data is crucial to obtain a more systematic representation of the state of a 
dynamical system in many physical applications. For most applications, observed 
data which should ideally be defined on the whole physical domain, can be mea¬ 
sured only discretely, often with inadequate resolution and contaminated by errors. 
The problem of constructing approximate solution of dynamic systems by fusing 
incomplete observed data with a mathematical model has been the focus of many 
researchers in the past decades. For example, designing feedback control algorithms 
for linear time-invariant systems has been the focus of Luenberger and many others 
since the early 60’s and 70’s (see e.g., [28] and references therein). The problem of 
constructing approximate solution for classes of non-linear-nth-order systems was 
later studied by Thau in [42]. A recent review by Nijmeijer [34] addressed the prob¬ 
lem of the synchronization from a non-linear control theory prospective. Analyzing 
the validity and success of a data assimilation algorithm when some state variable 
observations are not available as an input on the numerical forecast model, in the 
context of meteorology and atmospheric physics, was also studied in [8, 20, 21, 25]. 

The classical method that was proposed, see, e.g., [8] and [11], is to insert ob¬ 
servational measurements directly into a model as the latter is being integrated in 
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time. One way to exploit this is to insert low Fourier mode observables into the 
equation for the evolution of the high Fourier modes. This was the approach taken 
for the 2D Navier-Stokes (NSE) in [6, 24, 29, 36, 37]. The same approach was also 
taken by the authors in [23] to present a discrete-in-time data assimilation algo¬ 
rithm for the Lorentz and the 2D NSE systems. The authors in [23] and [36] were 
the hrst to study the data assimilation problem from a PDE prospective following 
ideas from the theory of determining modes for infinite-dimensional dynamical sys¬ 
tems. The previously mentioned theoretical work assumed that the observational 
measurements are error free. We refer to some concurrent related works that al¬ 
low for noisy observations in [5, 30] which were motivated by the earlier analytical 
studies in [23] and [36]. 

In [3], the authors introduced a simple finite-dimensional feedback control al¬ 
gorithm for stabilizing solutions of infinite-dimensional dissipative evolution equa¬ 
tions, such as reaction-diffusion systems, the Navier-Stokes equations and the Kuramoto- 
Sivashinsky equation. A computational study of this simple finite - dimensional 
feedback control algorithm was presented in [31]. Based on this control algorithm, 
a new continuous data assimilation algorithm was developed in [4]. In this new 
algorithm the coarse mesh observational measurements are incorporated into the 
equations in the form of a linear feedback control term. That is, the algorithm was 
built by adding to the model equations a term that nudges the solution towards 
the observations. In the 1970s (see e.g., [25] and references therein) this technique 
was called Newtonian nudging or dynamic relaxation method, analyzed typically 
in much simpler scenarios. The advantage is that, since no derivatives are required 
of the coarse grain observables, this works for a general class of interpolant opera¬ 
tors. The algorithm introduced in [4] was designed to work for general dissipative 
dynamical systems. 

The main idea can be outlined in terms of a general dissipative evolutionary 
equation 

| = ( 1 . 1 ) 

where the initial data u(0) = vP is missing. The algorithm is of the form 

dv 

— = F{v) - fi{Ih{v) - Ih{u)), (1.2a) 

u(0)=u°, (1.2b) 

where /x > 0 is a relaxation (nudging) parameter and is taken to be arbitrary 
initial data. Ih(') represents an interpolant operator based on the observational 
measurements of a system at a coarse spatial resolution of size h, for t G [OiTj. 
Notice that if system (1.2) is globally well-posed and Ih{v) converge to Ih{u) in 
time, then we recover the reference u(t,x) from the approximate solution v{t,x). 

The main task is to find estimates on /x > 0 and h > 0 such that the approximate 
solution v(t) is with increasing accuracy to the reference solution u{t) as more 
continuous data in time is supplied. 

It is worth noting that the continuous data assimilation in the context of the 
incompressible 2D NSE studied in [4] under the assumption that the data is noise 
free. The case when the observational data contains stochastic noise is treated in 
[7]. The authors in [7] established resolution conditions which guarantee that the 
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limit supremum, as the time tends to infinity, of the expected value of the L^- 
norm of the difference between the approximating solution and the actual reference 
solution, (i.e. the error), is bounded by an estimate involving the variance of the 
noise in the measurements and the spatial resolution of the collected data, h. 

It is also worthwhile to mention that the authors in [33] analyzed an algorithm 
for continuous data assimilation for 3D Brinkman-Forchheimer-extended Darcy (3D 
BFeD) model of a porous medium, a model equation when the velocity is too large 
for classical Darcy’s law to be valid. A similar algorithm for stochastically noisy 
data is at hand combining ideas from the work in [7] and [16]. Furthermore, the pro¬ 
posed data assimilation algorithm can also be applied to several three-dimensional 
subgrid scale turbulence models. In [1], it was shown that approximate solutions 
constructed using observations on all three components of the unfiltered velocity 
held converge in time to the reference solution of the 3D NS-a model. An abridged 
data assimilation algorithm for the 3D Leray-a model, using observations in only 
any two components and without any measurements on the third component of the 
velocity field, is analyzed in [17]. 

A particular application that we have in mind in this paper is hows in porous 
media, subjected to heating from below and cooling from the top. Flows in porous 
media are connected to many important problems in geophysics (see, e.g., [39] and 
references therein) as well as in biological systems and biotechnology (see, e.g., 
[43]). 

We consider D = [0, L] x [0, 1] x [0,1] to be a box in of porous media satu¬ 
rated with a huid. The side walls of the box are insulated, and the box is heated 
from below with constant temperature Tq and cooled from above with constant 
temperature Ti < Tq. After some change of variables and proper scaling, (see, e.g. 
[10, 39]), the governing non-dimensional equations of the convecting huid through 
a porous medium are: 

Qu 

j—+ u + Vp = Ra9k, (1.3a) 

— -Ae+(u-V)e-u-k = 0, (1.3b) 

at 

V ■ u = 0, (1.3c) 

u{0;x,y,z) =u°{x,y,z), e{0;x,y, z) = e°{x,y, z), (1.3d) 


subject to the boundary conditions: 

9{t;x,y,0) = 9{t;x,y,l) = 0, 

99, ^ , 89, ^ ,89, 89, , , „ 

'^y9,0,y,z) = —{t-,L,y,z) = —{t-,x,0,z) = —[t-x,l,z) = 0, 

and 


(1.3e) 

(1.3f) 


u • h = 0 , on 9D. (l-3g) 

Here, u{t]x,y,z) = {ui(t; x,y, z),U 2 {t; x,y, z),U 3 {t]x,y, z)) is the huid velocity, 
p = p(t; X, y, z) is the pressure. 9 = 9{t] x, y, z) is the scaled huctuation of the tem¬ 
perature around the steady state background temperature prohle (Ti — To)z + Tq 
and it is given hy 9 = T — (^r^^ — z), where T = T{t-, x, y, z) is the temperature of 
the huid inside the box D, Tq is the temperature of the huid at the bottom and Ti is 
the temperature of the huid at the top. The non-dimensional parameter 7 is equal 
to where Pr is the Darcy-Prandtl number representing a measure of the ratio 
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of the viscosity to the thermal diffusion coefficient. The non-dimensional parameter 
Ra is the Rayleigh-Darcy number, which is a measure of the ratio of the driving 
force (coming from the imposed average temperature gradient), to the viscosity and 
thermal diffusion in the system. The vector k is the unit vector in the z-direction 
and n = n(x, y, z) is the normal vector to the boundary dfl at the point (x, y, z). We 
remark that the temperature fluctuation 9{t; x, y, z) satisfies an advection-diffusion 
equation (1.3b). Also, the evolution of the fluid velocity u(t;x,y,z) and pressure 
p{t;x,y,,z) is described by Darcy’s law (1.3a), this replaces the conservation of 
momentum equations in the Navier-Stokes equations. 

The Benard convection in porous media problem (1.3) with 7 > 0 (namely, 
the Darcy-Prandtl number Pr is finite), was studied in [12] and [13]. There, the 
existence and uniqueness of global weak and strong solutions were established. It 
was also shown that the temperature 0(t] x, y, z) satisfies the maximum principle 
and has an absorbing ball in L^’(D) for each finite p. The authors of [14] showed 
that for initial data (it°, 0 °) € L^(D) x L°“(D), the system has a finite-dimensional 
global attractor A C L^(D) x L°°(D) which attracts in the x L^(D) metric. 

Moreover, they showed that the system has exponential attractors. The Gevrey 
regularity (spatial analyticity) of the solutions on the global attractor was later 
studied in [35]. There, a rigorous lower bound estimate on the radius of analyticity 
was obtained. 

When 7 = 0 (namely, the Darcy-Prandtl number Pr is infinite), system (1.3) 
was studied in [32] . All the constraints on the boundary are similar to the nonzero 
Darcy-Prandtl number case, while we note that in the 7 = 0 case, there is no need 
to specify the initial data m°(x, y, z). One can recover tt°(x, y, z) by solving equation 
(1.3a) given 9°{x,y,z). This is unlike the non-zero Darcy-Prandtl case where one 
has to specify both initial data u°{x,y, z) and 9^{x,y, z). It was shown in [32] that 
the system for the 7 = 0 case has global real analytic solutions and admits a real 
analytic global attractor. Moreover, it was also shown that the standard Galerkin 
solution of the system, when 7 = 0, converges exponentially fast, in the wave 
number, to the exact solution of the system. This in turn justifies the computational 
results reported in [ 22 ]. 

We propose a continuous data assimilation algorithm for the 3D Benard convec¬ 
tion in porous media employing coarse mesh measurements of only the temperature 
field. We incorporate the observables as a feedback term in the original evolution 
equation for the system. This algorithm can be implemented with a variety of 
finitely many observables: for instance low Fourier modes, nodal values, finite vol¬ 
ume averages, or finite elements. Our algorithm is given by the following system: 


+v + yq = Rayk, (1.4a) 

at 

dp 

— - A 77 -I- (u • V)?? -V ■k = -y{Ih{p) - lh{0)), (l-4b) 

V-?; = 0, (1.4c) 

w(0;x,y,z) = ?;°(x,y,z), p(Q\x,y,z) = rf{x,y,z), (1.4d) 
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subject to the boundary conditions: 



r]{t;x,y,0) = r]{t;x,y, 1 ) = 0 , 

(1.4e) 

|^(<; 0 ,y, 2 :) 


(1.4f) 


V ■ n = 0, on i9D. 

(l-4g) 


Here, g is a modified pressure, and we note that the initial conditions vq and 779 , may 
chosen arbitrarily, (for e.g., both initial conditions are set to zero). In this paper, 
we will only consider interpolant observables that are given by a linear interpolant 
operator —>■ satisfying the approximation property 

ll</5-4(<p)|lL=(n) < co/i||75||jj-i(o), (1-5) 

for every (p G where cq > 0 is a dimensionless constant. One example of an 

interpolant observable that satisfies (1.5) is the orthogonal projection onto the low 
Fourier modes with wave numbers less than 1/h. A more physical example are the 
volume elements that were studied in [27]. The algorithm and our analysis in this 
paper will carry on to the second type of interpolants treated in [4]. This second 
type is given by a linear operator Ih : —?► T^( 0 ), together with 

Il7> - 4((/j)||i2(j.j) < cih ||75||jji(n) + C 2 h^ , (1.6) 

for every (p G where ci,C 2 > 0 are dimensionless constants. An example of 

this type of interpolant observables that satisfies ( 1 . 6 ) is given by the measurements 
at a discrete set of nodal points in fl. We are not treating the second type of 
interpolants in this paper only for the simplicity of presentation. For full details on 
the analysis for the second type of interpolants we refer to [15] and [16]. 

Our current analytical investigation is motivated by some recent studies in [15] 
and [16]. In [15], a continuous data assimilation scheme for the two-dimensional 
incompressible Benard convection problem was introduced. The data assimilation 
algorithm in [15] constructs the approximate solutions for the velocity u and tem¬ 
perature fluctuations 6 using only the observational data, 4 (u), of the velocity 
field and without any measurements for the temperature (or density) fluctuations. 
Inspired by the recent algorithms proposed in [4, 15], we introduced an abridged 
dynamic continuous data assimilation for the 2D Navier-Stokes equations (NSE) in 
[16]. The proposed algorithm in [4] for the 2D NSE requires measurements for the 
two components of the velocity vector field. On the other hand, in [16], we establish 
convergence results for an improved algorithm where the observational data needed 
to be measured and inserted into the model equation is reduced or subsampled. 
Our algorithm there requires observational measurements of only one component 
of the velocity vector field. Ideally, one would like to design an algorithm, for the 
Benard convection model, based on temperature measurements only, but the anal¬ 
ysis to support or dispute the existence of such an algorithm is a subject of ongoing 
research. In fact, the authors of this paper believe that knowing the temperature 
trajectory does not, in general, determine the velocity vector field in the 2D Benard 
convection model (see [2] for a computational study supporting this claim). It is 
worthwhile to mention that earlier conjecture of Charney et. al. [ 8 ], that complete 
temperature history of the atmosphere for certain simple atmospheric models, will 
determine other state variables, was tested in [20, 21] (see also [2]). It was noted 
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in [21] that in practice, the complexities in weather forecast models together with 
numerous problems in the real data make it uncertain that assimilation of tem¬ 
perature data alone will yield initial states of arbitrary accuracy. The authors in 
[21] developed and tested different time-continuous data assimilation methods us¬ 
ing temperature data. Their results suggested that temperature data can have an 
impact on the numerical weather forecasts that is sensitive to the quality of such 
observed data and is highly affected by the assimilation method used. The authors 
concluded that designing more refined methods for data assimilation and develop¬ 
ing better numerical models can lead a major improvement in numerical weather 
prediction using temperature measurements alone. 

Our main contribution in the current study is to give a rigorous justification 
that a data assimilation algorithm based on temperature measurements alone can 
be designed for the 3D Benard convection in porous media. We provide explicit 
estimates on the relaxation (nudging) parameter p and the spatial resolution h of 
the observational measurements, in terms of physical parameters, that are needed 
in order for the proposed downscaling algorithm to recover the reference solution 
under the assumption that the supplied data are error free. In the case where the 
observational measurements are not error free, one can estimate the error between 
the solution of the algorithm and the exact reference solution of the system in terms 
of the error in the measurements. While the typical scenario in data assimilation is 
to choose /i depending on h, in our convergence analysis we choose our parameters 
^ and h to depend on physical parameters. More explicitly, we choose fj, to depend 
on the bounds of the solution on the global attractor of the system and then choose 
h to depend on The philosophy here is that in order to prove the convergence 
theorems, we need to have a complete resolution of the flow, so h has to depend on 
the physical parameters, such as the Rayleigh number (i?a) in this case. 

Numerical simulations in [18] and [19] (see also [23]) have shown that, in the 
absence of measurements errors, the continuous data assimilation algorithm for the 
2D Navier-Stokes equations performs much better than analytical estimates in [4] 
would suggest. This was also noted in a different context in [36] and [37]. It is likely 
that the data assimilation algorithm studied in this paper will also perform much 
better than suggested by the analytical results, i.e. under more relaxed conditions 
than those assumed in the rigorous estimates. This is a subject of future work. 

The outline of the paper is as follows. In section 2 we give some preliminaries. In 
section 3.1, we present the convergence analysis of our proposed data assimilation 
algorithm in the case of zero Darcy-Prandtl case and the nonzero Darcy-Prandtl 
case in section 3.2. 


2. Preliminaries 

For the sake of completeness, this section presents some preliminary material 
and notation commonly used in the mathematical study of hydrodynamics models, 
in particular in the study of the Navier-Stokes equations (NSE) and the Euler 
equations. For more detailed discussion on these topics, we refer the reader to, e.g., 
[9], [38], [40] and [41]. 

Let LP{Q) := iP : 17 —R and iJ^(I7) := : 17 —R be denote the usual 

Lebesgue space and H^-Soholev space, respectively, for 1 < p < cx) and k gR. We 
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define the spaces: 

V := {w G : M • n = 0 on dQ and V • u = 0 in f2}, 

V := {6 G C°°(rf) : 9 satisfies the boundary conditions (1.3e) and (1.3f)}, 

H := the closure of V in the norm, 

V := the closure of V in the norm, 

H := the closure of V in the LF‘{^) norm, 

V := the closure of V in the norm. 


We define the inner products on H and H by 

3 


(m, ic) = / UiWi dxdydz, and {9,r])= / Oijdxdydz, 

respectively. We also define the inner product on V and V by 

3 . ^ r 

{{u,w)) = / ^3 UidjWi dxdydz, and {{6,r])) =''^^ / djOdjTj dxdydz 

■ ■ _T O •_ ^ 




respectively. Note that, thanks to the boundary conditions (1.3e) and (1.3f), H-Hy = 
((•c))^/^ = ll V-||j;^ 2 (q) is a norm on V due to the Poincare inequality (2.2), below. 
We define the Helmholtz-Leray projector as the orthogonal projection from 
onto H and define A = —A subject to the boundary condition (1.3g) and 
(1.3f) with the domain 

ViA) = {9 G : 9 satisfies (1.3e) and (1.3f)}. 


Using the Lax-Milgram Theorem and the elliptic regularity in the box fl, the linear 
operator A is self-adjoint and positive definite with compact inverse : H — 'D(A). 
Thus, there exists a complete orthonormal set of eigenfunctions Wi in H such that 
Awi = XiWi where 0 < < Ai+i for i G N. 

We define the bilinear from •) : V x 'D{A) H, such that 

B{u,9) := (m-V)6 l, 

for every it G V and 9 G PiA). Using the boundary conditions (1.3e) and (1.3f), 
one can easily check that 


{B{u,9),9)=0, 


( 2 . 1 ) 


for every it G V and 9 G P{A). 

We recall the Poincare inequality: 

for all ^gU, (2 .2) 

where Ai is the smallest eigenvalue of the operator A. 

Let y be a Banach space. We denote by LP[\fi,T\,Y) the space of (Bochner) 
measurable functions t i->- 'w{t), where w{t) G Y, for a.e. t G [0,T], such that the 
integral fj" ||w(OIIf is finite. 

Furthermore, inequality (1.5) implies that 


11 ^ ^h{9)\\j^2^Q^ 


< coh 


A^/^9 


(2.3) 


for every 9 gV . 


L^(n) ’ 
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We recall the following existence and uniqueness results for the 3D Benard con¬ 
vection problem in porous media (1-3). Hereafter, c will denote a universal dimen¬ 
sionless positive constant that may change from line to line. 

Proposition 2.1. [32] Given 6 G V, there exists a unique solution u G V of the 
problem (1.3a) and (1.3c), with 7 = 0 , subject to the boundary condition (1.3g). 
Moreover, u satisfies 

||w||y-^ ci?a ||0||y . (2.4) 

Theorem 2.2. [32] Let 9^ G V, then system (1.3), with 7 = 0, has a unique global 
strong solution {u, 9) that satisfies: 

uGC{[0,T]-,V)nL\[0,T]-,{H\n)f), and ^ g 

9GC{[0,T]-,V)nL^{[0,T]-V{A)), and ^ G L^{[0,T]-, H), 

for any T > 0, and depends continuously on the initial data 9^. Moreover, the 
solution satisfies the following bound: 

limsup|| 6 »(t)||^=o(s. 2 ) ^ (2.5) 

t—¥oo 

and the system admits a compact finite-dimensional global attractor .A C V x H. 

Theorem 2.3. [lj[ //m° G H and 0° G then for any T > 0, system (1.3), 

with 7 > 0, has a unique weak solution {u, 9) that satisfies 

uG L°°{[0,T]-U), and 9 G L°°{[0,T]-L°°{n)) D L^{[0,T]-V). 
Moreover, for every finite p, 

limsup||6l(t)||^P(j.j) < (2.6) 

More precisely, 6{t;x,y,z) can be decomposed such that 

S{t;x,y,z) = 9i{t-,x,y,z) + 92{t-,x,y,z), (2.7) 

where —1 < 9i{t', x,y, z) < 1, and that there exist two positive constants a\ and 02 
depending on p and LI such that 

ll^ 2 (t)|lLr(n) ^ “le 

If the initial data 9^ satisfies m < 9^{x,y,z) < M, a.e. in fl, for some m > —1 
and M < 1, we have 

m < 9{t-,x,y,z) < M, a.e. in LI, (2-8) 

for all t > 0. 

Finally, if € V and 9^ G L°°{Ll), then there exists a unique strong solution 
{u, 9) that satisfies 

uGC{[0,nV), and ^ g L^{[0,T]-,U), 

9GL°°{[0,T];L^{Ll))nL\[0,T];V{A)), and ^ g L\[0,T]; H), 

and the system admits a compact finite-dimensional exponential attractor A C Vx 

y. 
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3. Convergence Results 


In this section, we will establish the global existence, uniqueness and stability of 
solutions of system (1.4), when the observable data satisfy (1.5). 

Since we assume that (m, 0) is a reference solution of system (1.3), then it is 
enough to show the existence and uniqueness of the difference {w, = {u — v, O — rf). 

In the proofs of Theorem 3.1 and Theorem 3.2 below, we will drive formal a-priori 
bounds on the difference (w,^), under the conditions that p, is large enough and h 
is small enough such that pcgh < 1. These a-priori estimates, together with the 
global existence and uniqueness of the solution {u,6), form the key elements for 
showing the global existence of the solution (u, ij) of system (1.4). The convergence 
of the approximate solution (u, p) to the exact reference solution (u, 0) will also be 
established under the same conditions on the nudging parameter p stated in (3.1) 
and (3.9), when 7 = 0 and 7 > 0, respectively. Uniqueness can then be obtained 
using similar energy estimates. 

The estimates we provide in this section are formal, but can be justified by the 
Galerkin approximation procedure and then passing to the limit while using the 
relevant compactness theorems. We will omit the rigorous details of this standard 
procedure (see, e.g., [9, 12, 14, 32, 38, 41]) and provide only the formal a-priori 
estimates. 


3.1. The infinite Darcy-Prandtl number case (7 = 0 ). 

Theorem 3.1. Let satisfy the approximation property (1.5) and {u,6) be a 
strong solution in the global attractor of (1.3) with 7 = 0 . Suppose that G H 
and G H. Assume that p > 0 is large enough such that 

p + ^ > 2cRa^ + 4:Ra, (3.1) 

and that h > 0 is small enough such that pc^h^ < 1. Then, there exists a unique 
weak solution (v,p) of system (1.4) with 7 = 0 such that 

uGC([0,T];H)nL2([0,T];(V), and ^ g L\[0,T];V'), 

pGC{[0,T]-,H)nL^i[0,T]-,V), and ^ G L\[0 ,T]-v'). 

Moreover, ||u(t) — > 0, and \\0{t) — p{t )\\^]^2 — >■ 0, at an exponential 

rate, as t —>■ 00 . 


Proof. Define w := u — v and ^ := 0 — p. Then, in functional settings, w and f will 
satisfy the equations: 


w = RaPa{fk), 

^ + B{v,f) + B{w,9) -wk = -plhiO- 


(3.2a) 

(3.2b) 


Taking the 77-inner product of (3.2b) with we obtain 

2 


2^ll?ll^(0) 


^ 1 / 2 ^ 


L^(n) 


+ iBiv,ao + mw,0),o 

= {wk,^)- p{lh{0,0- 
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Thanks to (2.1), we have 

We also notice from (3.2a) that 

w = RaPai^k), in t 2([0,T];H), 
for any T > 0. This implies that 


(3.3) 


{w ■ k,^) = Ra{Pa{^k),^) < Ra Pa{^k) 
Using and Holder’s inequality, we get 

|(S(ii;, 0), 01 = I(S(«^, 0,0)1 

< \\e\ 


L^(n) 


IICIlL^(n) —ll'CllL^(n) ■ (3-4) 


lRq) iriiL“(n) 


yll/2^ 


Lqo) 


< 


yll/2^ 


L^(n) 


llL“(n) 


c||6»||roo^o^ l|lC||r2 


2 

L^(0) 


Thanks to the equation (3.2a) we have 

I(S(^,0),OI<^ 


L-’(n) 


cRa^ 


lL“(n) ll^llL=(n) 


(3.5) 


The approximation inequality (2.3) and Young’s inequality imply 

-m(4(0, 0 = -m( 4(0 -00“/^ ll^llL’'(n) 

< /I l|4(0 ~ 'fllL’’’(n) ll^llL’^(n) ~ M ll^llL’’’(n) 


< 


/iCg/l^ 


yll/2^ 




L^{n) ■ 


(3.6) 


Using assumption /rcg/i^ < 1 and estimates (3.3)-(3.6), we conclude that 


The Poincare inequality (2.2) yields 
d „.„2 A 


^ lieil^(n) + Y lieil^(o) < {cRa^ l|0fL<==(n) + 2i?a - n) ||OI^(a) ■ 

We define 


Then, 


Y{t) := fi+^- 2Ra - cRa^ ||0|li=°(n) ■ 


Yll^ll^(O) + «Wlieil^(O)<0. 


(3.7) 


The uniform bound (2.5) implies that: for any fixed £ > 0, there exists a time 
to(£) > 0 such that 

l|0(OllL°“(n) < 1 + £, 

for all t > to. Taking £ = 1, then there exists a time to > 0 such that ||0(t)||£°o(Q) < 
2, for all t > to- Then, we have 

a{t) > ^ + Y ~ 2i?a — cRa^, 
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and 

a{t) < fJ. + -^ + 2Ra + cRa^ < oo, 

for all t > to- Now, the assumption (3.1) implies that 

a{t) > 2Ra + cRa^ > 0, 

for all t > to- 

By Gronwall’s Lemma, it follows that 

Pit) - vit)\\l^(^n) = IICWIli=(n) ^ 0, 

at an exponential rate, as t —?> cx). Finally, the equation (3.2a) yields that 

II^WIlL=(n) < ^«IIC(i)lli2(o), 

thus, 

Mt) - u(t)||^2(j-j) = 11^11^2(5-2) ^ 0^ 

at an exponential rate, as t —>■ oo. 


(3.8) 


□ 


3.2. The finite Darcy-Prandtl number case (7 > 0). 

Theorem 3.2. Fix 7 > 0, and (u, 0) be a strong solution in the global attractor 
of (1.3). Suppose that G H, rj^ € H and that R satisfies the approximation 
property (1.5). Let p > 0 be large enough such that 

2/i + Ai >2^^+2c7(l + AC^)^ (3.9) 

7 

and h > 0 is small enough such that < 1. Then, there exists a unique weak 

solution iv,rj) of system (1.4) corresponding to the same 7 > 0 such that 

uGC([0,T];H), and ^ l\[0,T]-,v') 

r,€Ci[0,T]-H)nL^{[0,T]-,V), and ^ g L^([0,T];V'). 

Moreover, ||u(t) — v{t )\\\^2 + ||0(t) “ 0 (t)||^ 2 ( 52 ) —0, at an exponential rate, 
as t ^ 00 . 

Proof. Define w := u — v and ^ := 0 — p. Then, in functional settings, w and ^ will 
satisfy the equations: 

7 ^ + w = RaPai^k), 

+ 0 + B{w, e)-wk = -plhi^. 


(3.10a) 

(3.10b) 
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Taking the H-inner product of (3.10a) with w and the 77-inner product of(3.10b) 
with respectively, we obtain 

ll^llL=(n) + ll«^llL=(n) = RaiC,w- k), 


Thanks to (2.1), we have 
Also, estimate (3.6) shows that 


= {wk,0- 

mv,o,o = o- 

„l/2^ ^ M ll^l|2 


-nihiO.O < ^ f lieil^(n) ■ 


(3.11) 

(3.12) 


By Holder’s inequality, we have 

|(i3(w;, 0), 01 = I(S(«^, 0,0)1 

< IkllL^(n) l|0|lL==(n) 


L^(n) 


< 


1 


^1/2^ 


L^{U) 


c||0| 


L“(n) II^IIl^(O) 


(3.13) 


The above estimates (3.11)-(3.13), the Cauchy-Schwarz inequality, and the Poincare 
inequality (2.2) yield 

^ ^ = RaiX.w k) 


2 ^ ii“'iiLqn) "E ii‘"iiL^(n) 


- 2 + 2 ll^llL’^(n) 




2 


L\n) + 2 - 2 


^1/2^ 


L2(n) 4 

2 


1 / 2 
^1/2^ 


2 


+ c l|0|li~(n) Ik||i2(j2) + cAi ^ \\w\\l2^^^. 

Using the assumption < 1 imply 

^ IkllL’^(n) + - Iklli^cn) < — IICIli2(n) > 


(3.14a) 


^ ll^llL=(n) + (/^+ "T 


iL'2{a) ^ ^ + ||0|li~(f2)^ ll'^llL^(n) • (3.14b) 


Multiplying (3.14a) with ||'R'||^ 2 (-j. 2 ^ and (3.14b) with ||0lL2(n), respectively, we 
have 
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Using the Cauchy-Schwarz inequality, we get 



The above inequality can be rewritten as 


I (ikll^(n) + lieil^(o)) + (ikll^(n) + < 0, (3-15) 

where 

a{t) = minji, 2 /r + Ai - - cy |. 

Recall that, by Theorem (2.3), the solution 0{t; x, y, z) of the 3D Benard problem 
in porous media satisfies the Maximum Principle and 

limsup||6»(t)||^=o(s^) < 1. 

t—^OO 

Then, by a similar argument as in the proof of the previous theorem, we conclude 
that there exists a time to > 0 such that 


and 


a(t) > min ■( —, 2 /r + Ai — - 07(1 + A;^ 

.7 7 


y{t) < min <{ -, 2 ^ + Ai + + 07(1 + A;^ < 


for all t > tg. Now, assumption (3.9) implies that 


a{t) > min i + 07(1 + A;^ ^)^ f- > 0 , 


,7 7 

for all t > to- 

By Gronwall’s Lemma , it follows that 


Mt) - + ll^(i) - 0^ 


at an exponential rate, as t —>■ 00 . 


(3.16) 


□ 


Remark 3.3. Suppose that initial data is contained in double the size of 

the absorbing ball (for example, one {v°,r]°) = (0,0)). Consider the approximate 
solutions of system (1.4) with interpolant operators Ih satisfying either of the ap¬ 
proximate properties (1.5) or (1.6). Then one can find /i > 0 large enough and h > 0 
small enough, depending on the physical parameters of the underlying system, such 
that v{t) converges in the strong sense to the reference solutions of system (1.3), i.e. 
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||u(t) — L'(t)||v + 11^(0 “ ^(Ollv ^ O’9''^ exponential rate, as t —)► oo. The proof 
of this statement follows similar adjusted arguments and ideas to those presented 
in [15] and [16]. 
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